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Abstract 

Background: Delineation of the target volume is a time-consuming task in radiotherapy treatment planning, yet 
essential for a successful treatment of cancers such as prostate cancer. To facilitate the delineation procedure, the 
paper proposes an intuitive approach for 3D modeling of the prostate by slice-wise best fitting ellipses. 

Methods: The proposed estimate is initialized by the definition of a few control points in a new patient. The method 
is not restricted to particular image modalities but assumes a smooth shape with elliptic cross sections of the object. 
A training data set of 23 patients was used to calculate a prior shape model. The mean shape model was evaluated 
based on the manual contour of 1 0 test patients. The patient records of training and test data are based on axial 
T1 -weighted 3D fast-field echo (FFE) sequences. The manual contours were considered as the reference model. 
Volume overlap (Vo), accuracy (Ac) (both ratio, range 0-1, optimal value 1) and Hausdorff distance (HD) (mm, optimal 
value 0) were calculated as evaluation parameters. 

Results: The median and median absolute deviation (MAD) between manual delineation and deformed mean best 
fitting ellipses (MBFE) was Vo (0.9 ± 0.02), Ac (0.81 ± 0.03) and HD (4.05 ± 1 .3)mm and between manual delineation 
and best fitting ellipses (BFE) was Vo (0.96 ± 0.01 ), Ac (0.92 ± 0.01 ) and HD (1 .6 ± 0.27)mm. Additional results show a 
moderate improvement of the MBFE results after Monte Carlo Markov Chain (MCMC) method. 

Conclusions: The results emphasize the potential of the proposed method of modeling the prostate by best fitting 
ellipses. It shows the robustness and reproducibility of the model. A small sample test on 8 patients suggest possible 
time saving using the model. 

Keywords: Delineation, Ellipse model, Empirical Bayes, Prostate, Radiotherapy treatment planning, Statistical shape 
analysis 



Background 

Prostate cancer is the second most diagnosed cancer 
accounting for 14 percent of all cancers diagnosed world- 
wide [1]. It is most common in males over the age of 50, 
and has the highest incidence rate in the developed coun- 
tries. Aggressive tumors are usually treated with extern 
radiotherapy or brachytherapy which requires a precise 
treatment plan for the target volume. In any type of radio- 
therapy treatment, radiation of healthy tissue should be 
minimized while maintaining the desired dose to the tar- 
get volume. Therefore, a successful treatment of prostate 
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cancer relies on an accurate segmentation of the prostate 
from the surrounding tissue, by image-based description 
of the shape and location of the target volume. The vol- 
ume of interest is characterized by a smooth shape, and 
for this reason an algorithmic description of the volume is 
feasible. 

Transrectal ultrasound (TRUS), magnetic resonance 
(MR) and computed tomography (CT) images are the 
three main imaging techniques used in diagnosis, treat- 
ment planning and follow-up examination of prostate 
cancer. Smith et al. [2] investigated the effects of these 
imaging techniques on the properties of the prostate 
volume. A collection of methods available for prostate 
segmentation is reviewed by Ghose et al. [3]. In addi- 
tion to the methods presented by Ghose et al, alternative 
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approaches are available in the literature, such as the 
medial or skeleton representation of the prostate [4-8]. 
The present work proposes a segmentation method which 
falls into the category of deformable meshes in Ghose et al 
[3], but refers to the term geometrical parametrization as 
described in Dryden and Mardia [9]. The main focus of 
this paper is the development of a statistical shape model 
for the prostate. An overview about this type of models in 
3D medical image segmentation is presented for example 
by Davies et al. [10] and Heimann and Meinzer [11]. 

The works of Saroul et al. [12] and Mahdavi et al. [13] 
are related to the stacked ellipses parametrization method 
used in this paper. Mahdavi et al. [13] proposes a 3D 
ellipsoid shape of the prostate in warped transrectal 3D 
ultrasound images based on control points. This method 
extends the warping idea proposed in Badiei et al. [14] 
from 2D to 3D ultrasound images. On the contrary, we 
focus on slice-wise best fitting ellipses which will intro- 
duce more flexibility into the model, e.g., between the 
positions and lengths of the first and second axes of the 
ellipses between neighbor slices. The approach of slice- 
wise best fitting ellipses has similarities to a tubular medial 
representation [15]. 

Beside the single segmentation of the prostate, several 
attempts have been tried out for a joint segmentation of 
neighbor organ and structure to gain improved segmenta- 
tion results [16-18]. 

To our knowledge, despite the substantial effort in this 
area, no widely implemented algorithm exists. In oncology 
departments this means that the physician has to delin- 
eate the prostate slice by slice. This is time-consuming 
and inefficient. We propose a less ambitious approach 
compared to more sophisticated models, such as skele- 
tal models as discussed above, in that we use a method 
that gives a useful starting point for the physician after 
the definition of few control points. Given the initial esti- 
mate of the volume of interest, the physician can adjust the 
estimate according to their evaluation of the image rather 
than starting from scratch. By this approach, we obtain 
the same accuracy with less effort. The main points in our 
approach are as follows: First, we accept that the algorithm 
cannot give a fully precise description of the volume. Our 
main aim is therefore to give a good estimate which can 
be used as a starting point for the physician. Second, we 
use a simple ellipse model that is easy to interpret and 
understand. Our hypothesis is that a more efficient use 
of physicians in Radiotherapy Treatment Planning (RTP) 
of patients with prostate cancer can be obtained by an 
easy-to-interpret semiautomatic tool. 

Figure 1 shows an example of the initial estimate we 
typically obtain for a single image slice. The dashed line 
in (a) to (e) describes the manual contour while the solid 
line shows the best fitting ellipse including the two prin- 
cipal axes for the observed data of this slice. Note that the 



fitted model is very much in agreement with the manual 
line, indicating that the stacked ellipses model gives a good 
description of the object of interest. The solid lines in 
(f )-(j) shows the outcome from our model in this situation 
together with few defined control points. 

This result shows a typical performance of the method, 
and that the estimate is close to the best fit we can obtain 
with the ellipse model. The full processing demands little 
computational resources, such that the suggested delin- 
eation can be presented immediately. The example is dis- 
cussed further in the Methods and Results and discussion 
section. 

The rest of the paper is organized as follows. In the 
Methods section, we introduce the data sources and the 
proposed stacked ellipses model, and discuss the shape 
space and statistics along with constraints and parame- 
ters. Results are presented in the Results and discussion 
section using a test data set to show the potential of the 
mean shape model, followed by a Conclusions section. 
Additional file 1 with further detailed discussion is avail- 
able online. 

Methods 

Preliminaries 

Each prostate must be described by a shape model in 
order to calculate statistics, e.g., by stacked ellipses as a 
parametric shape model. The parameters of a paramet- 
ric shape model can be estimated from a training set. 
The training set models also the geometric variability of 
anatomical structures by a shape probability distribution. 
The training set contains volume and contour informa- 
tion of segmented prostates from N patients. The volume 
information describes the image modalities (e.g., CT or 
MR) and the contour information the volume of interest 
as defined in the following. 

The volume information of each training set n = 1, 
. . . , N is defined by a 3-dimensional matrix V n where 
V n (i,h) contains the observed gray level in voxel (/,/?), 
i = i 2 ) e {1, . . . ,h} x {1, ... ,72} are the pixel indices 
in a slice, where typically h = h> and h e {1, ... ,77} is 
the number of slices per data set. The number of slices 77 
is not necessarily the same for all patients in the training 
data sets. Therefore, we indicate 77 by H n and in the same 
manner 7i by I n \ and I2 by I n 2, but for simplicity we use 77, 
7i and I2 if the meaning is clear. 

In addition to the volume information, each training 
set n = 1, ...,N consists of contour information of 
the prostate, manually drawn by a physician. The con- 
tour information can be modeled by a (M x K n ) con- 
figuration matrix X n := (X n \, . . . ,X n K n ) with X n k = 
( x ik> x 2k' x 3k) T G R3 ' /c = !>..., K n , where K n defines 
the total number of available contour information points 
in a data set and M = 3 defines the dimension. We 
assume the contour information for an object is defined 
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(a) slice 6 



(f) slice 6 



(b) slice 9 



(c) slice 12 



(d) slice 15 



(g) slice 9 



(h) slice 12 



(i) slice 15 



(e) slice 17 




(j) slice 17 



Figure 1 Selected slices of MR data set 3 from the test data set. (a)-(e) Manual delineation of the prostate (dashed line) and best fitting ellipse 
(solid line), (f)-(j) Manual delineation of the prostate (dashed line), deformed mean shape (solid line) and defined control points in the first, center 
and last slice. 



in a sequentially sorted number L n of equidistant slices 
whereas each contour slice contains K n \ contour points, 
/ = 1, ...,L n . Hence it follows K n = J2i^nl an d 
X n — (X n i, . . . ,X n L n ). The image information in slice / is 
denoted by S n i and S n = {S n i, . . .,S nLn \ C V n and X n i 
defines the configuration matrix in slice S n i. 

In summary, the training population is given by the set 
{V, X}, with a set of volume information V = { V\, . . . , Vn} 
and configuration matrices X = {X±, . . . ,X/v}. We assume 
X n defines the configuration matrix for the correspond- 
ing data set V n and matches the volume information V n 
exactly. 

The contour information is often defined in a Patient 
based Coordinate System (PCS) whereas the volume 



information is defined in an Image based Coordinate 
System (ICS). The ICS can be transformed to PCS 
by a transformation matrix Adcm> which transform an 
image coordinate p im = (ii,i2>h) T to patient coordinate 
pP — (x,y,z) J \ The definition of Adcm and the relation 
between PCS and ICS (see Figure 2) is discussed in detail 
in the Additional file 1. In addition, we introduce a de- 
rotated PCS where volume and contour information are 
aligned to each other. 

Modeling 

The prior information inferred from the training set is 
incorporated into a shape model. We assume a stacked 
ellipse model as a shape prior for the prostate. Specifically, 



PCS 




de-rotated PCS 




ICS 


e.g. manual contour 




e.g. {dBFE nll dCP n } 




e.g. volume data 




An 

IV dCP 










sample space 








e.g. 


{BFE n i,CP n }, mean shape 





Figure 2 Visualization of different coordinate systems with example data. PCS: Patient based coordinate system (manual delineation line). 
ICS: Image coordinate system (volume data), de-rotated PCS: de-rotated patient based coordinate system with same scale and origin as PCS but 
same orientation as ICS (de-rotated best fitting ellipse dBFE n i, de-rotated control points dCP n ). sample space: The transformation matrix A n dcp maps 
the de-rotated data {dBFE n ,dCP n } to {BFE n ,CP n } in the sample space. 
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the prostate outline in slice S n i, I = 1, . . . , L n , n = 1, . . . , N 
is modeled by a slicewise best-fitting ellipse, as visual- 
ized in Figure 3. An ellipse in slice S n i can be uniquely 



described by p nl = (6 nl , a nL , (j) nL ) 
with 



l\ x (- 



7T_ 7T_ 1 

' 2> 2 J 



nl 



^ 2 defines the center in 



• position 0 

slice S n i, 

• length of principle axes a nl = (af, ot^) 1 G R% and 

• rotation angle (j) nl G (—§>§] in slice S n [. 

The rotation parameter 4> nl is defined corresponding to 
the ICS with origin 0 nl in slice S n [. The boundary of an 



ellipse p n i centered at 0 



nl 



in slice S n i is defined by 



C(Pnl) = \Rx- 



-\nl , 



(af) 2 



= 1,X G ' 



(1) 



and R = 



cos<p nl — sin0 



nl 



sin (j) nl cos (j) 



nl 



is a rotation matrix in M? with rotation angle (f) nl and # = 
(*i,* 2 ) r . 

The shape model described in this section requires the 
best fit of an ellipse C(p n i) to the contour information X n [ 
in each slice, i.e., we model X n \ = C(p n i) + e where e is an 
error with mean zero. The best-fitting ellipses provide us 
with a slice-by-slice parametrization of the prostate for all 
slices in each training shape. 

The problem of fitting an ellipse to geometric features 
like the contour is discussed widely in the literature (e.g., 
[19,20]). This work follows Ahn et al. [19], who proposed a 
least-square minimizer for X n [. The nonlinear estimate of 
parameters p n i 
minimize the error 



= (6»f,^ / ,off,^ / ,0 w/ ) r given X nl must 



g(Pnl) = (^nl - C(Pnl)) {^nl C(Pnl)) 



where C(p n i) is a set of nearest orthogonal points of X n i to 

C(Pni). 

Definition 1 (Best fitting ellipse (BFE)). A best fit- 
ting ellipse for slice S n i is defined by the set BFE n i := 
(O n Ka n K(i) nl ) T G R 2 x R2 x (-§,§], I = l,...,L m 
n — 1, . . .,N and minimizes the error function g, i.e., 
BFE nt = p nt with 



g{p n l) = min _ g(Pnl)- 

Pnl eR 2 



2 + x(" 



(2) 



The first and second principal axes must be reordered 
after calculation oiBFE n — {BFE n \ y . . .,BFE n i n \ in order 
to establish correspondence between parameters of adja- 
cent slices and across the population. Improved corre- 
spondence will support accurate statistics. The basic idea 
in our reordering procedure is to carry out the reordering 
corresponding to the lowest rotation angle of both prin- 
cipal axes to the first principal axis of the neighbor slice 
where the center slice is chosen as the basis. The rota- 
tion between the center slice M and an arbitrary slice is 
constrained by max(|0* — (p M \) = 7T,i G {1, ... ,1} after 
reordering. Therefore, the set BFE n of reordered best- 
fitting ellipses is an element of (R 2 x x (— tt, n]) Ln . 

A further improvement of correspondence is achieved 
by the introduction of two additional constraints in the 
parameter model. 

First, we relax the rotation parameter (p nl in case of cir- 
cularity. If both principal axes have the same length, the 
orientation of an ellipse is undefined. Therefore we penal- 
ize (p nl in the case of high circularity by taking § nl ' from 
the neighboring slices into account. Second, smoothing is 
performed between neighboring slices to avoid large for- 
ward and backwards rotations between 4> n ^ l ~ l \ cj) nl and 
-phe reordering algorithm and implementation of 
constraints are described in detail in Additional file 1. 
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The current implementation assumes the definition of 
control points CP n in the training data set{V n , X n , BFE n }, 
where BFE n e (R 2 x IR 2 x (— 7T,7t] is a reordered set 
of best fitting ellipses, n = 1, JV. Furthermore, the con- 
trol points have to be defined manually by a physician 
in a new patient data set. The control points are used 
to make the best fitting ellipses BFE n comparable and 
to transform the parametrized ellipses model to a com- 
mon position, scale and orientation by a transformation 
matrix A^ cp . The transformation matrix A^ cp maps the 
de-rotated prior data {dBFE n , dCP n ] to {BFE m CP n ] in the 
sample space, as depicted in Figure 2. In this article, we 
assume 6 control points in the first, center and last slice at 
the boundary of the prostate, i.e., 

CP n = {A n l ,...,Al,Pl,...,Pl,B n l ,...,Bl) 

as visualized in Figure 3. In addition, we have tested alter- 
native control point configurations. They are described 
together with the construction of A^ cp in Additional 
file 1. 

After transformation we have obtained a reordered and 
comparable set of best fitting ellipses 

BFE n = {BFE n i, . . . , BFE n L n } 

W\thBFE n i = (6 n \a nl ,4) nl ) T ,n = l,...,N,l= l,...,L n . 
The statistical analysis of the training data requires an 
equal number L\ = . . . = Lm to establish correspondence 
between the parameters of the best fitting ellipses. There- 
fore, we interpolate the set BFE n i to a common number 
L. 

When L is chosen, interpolation is done by indepen- 
dent cubic interpolation in each dimension, i.e., we find 
points of a one-dimensional function that underlies the 
data 0i ,0 O ,6o and 0 . The final interpolated 

best fitting ellipses are denoted by 

iBFE n = {iBFE n i, . . . , iBFE nL } . (3) 

These ellipses are used for the statistical analysis and 
computation of a mean shape model. To keep things sim- 
ple, we denote such a reordered, transformed and interpo- 
lated set of best-fitting ellipses by BFE ni = (0 nl , a nl , (j) nl ) T 
for the number L of contour slices with / = 1, ...,L 
and n = 1,...,N. The comparable set of best fit- 
ting ellipses BFE n is an element of the shape space 
(M 2 x M 2 x (-7r,7r]) L . 

Statistical analyses 

After reconstruction of our shape space we estimate the 
expectation and variance of the parameters of a mean 
shape model /jlbfe = {/^bfe 1 • • • » ^bff) w * tn I^bfe = 
(fig, /jL l a , from the training set BFE n i, I = We 



denote the mean shape mean best fitting ellipses (MBFE). 
In addition to the described ellipse parameters we define 
the position 6 nl = (0f,0% l ,0% l ) T in terms of a distance 
vector rj nl of 6 nl to a center curve defined by the control 
points. We model 6 l = £ / + rj l , where % l is analytically 
defined by L intersection points of the curve within each 
slice. Thereby, we are describing the mean shape which 
is closest to the control points. This approach is reason- 
able under the assumption that the control points are 
well defined. In Additional file 1 we explore various ways 
of describing the position parameter for different control 
point methods. 
The mean curve of the expected location is given by 

1 N 

^j = xY,°j l > U,2,3}, (4) 

i—l 

where \x l Q = (fig , a4 2 > ^ 3 ) r , I = 1,...,L. The variance 
and covariance are estimated by 

1 N 

K )2 = ]^ " ^ 7 ) 2 > {1,2, 3}, and (5) 

i—l 

i N 

*'o = N—[Y,Q U -riW U -ri) T - (6) 

i—l 

The length parameter is modeled by a log-normal distri- 
bution because a e . Thus we estimate the mean and 
variance of a = log(a) e IR 2 . The estimation of means and 
variances of the remaining parameters a, 0, r\ is according 
to (4-5). 

Following Dryden and Mardia [9] we suggest a prior 
distribution for a new data set as 

4~Af(/4>(<) 2 )^ 

a[ ~ log-TV (<) 2 ) with a[ = log(c^), 

4~Af(/4.(<) 2 )«=> 

a l 2 ~ \og-N ( M i 2 , (<rj 2 ) 2 ) with a l 2 = \og(a l 2 ), 
4> l -N^Aa^ 2 ), 

I — 1, . . . >L. If 0 l is defined according to the center curve 
given by the control points as described above, we model 

r]\ ~ N ^/x^., (tf^) 2 ) > * = 1» 2. Since the rotational param- 
eter is expected to have small variance it is not necessary 
to apply a circular distribution, and we assume normality. 

After constructing the shape model we estimate the best 
fitting ellipse BFEi parametrized by pi = (6 l , a 1 , (p l ) T , I = 
L in a new data set given the control points CP. This 
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is obtained through the posterior n(p | S) where su g5c 
V is the volume information and i = e I(p) is a 

set of indices within the ellipses p. The control points CP 
are used to deform the prior model tt(p). Therefore we 
model the posterior by an empirical Bayes approach [21]. 
The posterior 



7T(p I S, CP) oc L(S | p) * 7T(p I CP) 



(7) 



defines the posterior density of the deformed template 
7i (p | CP) given the the observed image. The Likeli- 
hood or image model L(S \ p) is the joint probability 
density function of the gray levels given the parametrized 
object p\cp> while p\cp defines the ellipses p deformed 
by the control points CP. The prior tt(p | CP) mod- 
els realistic variations from our mean shape \ibfe £ 
(IR 2 x IR 2 x (—tt,7t]) l given the control points. We are 
estimating the posterior distribution using a Markov chain 
Monte Carlo (MCMC) approach. The method and results 
are discussed in detail in Additional file 1. 

Evaluation 

We have evaluated the proposed method using 33 patient 
case studies. The training data set consists of N = 23 Tl- 
weighted Fast Field Echo (FFE) 3D Magnetic Resonance 
(MR) data. The mean shape model and variance is calcu- 
lated from the training data set and applied to a test data 
set of 10 MR FFE case studies. The splitting in test and 
training data is done according to the sequence of data 
acquisition. Each data set consists of H n Digital Imaging 
and Communications in Medicine (DICOM) image files 
and one DICOM region structure file, while the contour 
information of the prostate is stored in the header of a 
DICOM file without any image information. The voxel 
size (l x ,ly,lz) is (0.559mm, 0.559mm, 3mm) with a slice 
distance of 3.3mm of the data sets V n . Each slice con- 
sists of 288 x 288 voxels. The average number of slices 
with manual prostate contour information is 10.478 =b 
2.626 (meanzbstandard deviation) in the training set and 
10.5 =b 2.799 in the test set where the mean is given 
by \ii = ^ Yln=l Ln and the standard deviation by 

( W^l Yln=l (Ln — I^l) 2 ) ^ ♦ Figure 1 illustrates test patient 3, 
whose image set consists of 24 MR FFE slices whereas 12 
slices contain contour information. Slice 6 is the first slice 
where contour information of the prostate is available and 
the last slice is 17. 

An ethics approval was not required for this study 
under Norwegian law because the aim of this study was 
to develop a tool and not to obtain new knowledge on 
medicine or diseases. The study uses solely data that were 
collected during routine medical treatment independent 
of this study at the University hospital Northern Norway 
(UNN). Data was provided after a full anonymization 
following the required guidelines at UNN. 



Three metrics are used to compare the manual and the 
semi-automatic contours. In the axial slices, where the 
expert manual delineations are present, we calculate the 
Hausdorff distance (HD) by 

d^iX, Y) = max | max mind (#,;/), max min d(x, y) \ . 

I xeX yeY yeY xeX J 

(8) 

The Hausdorff distance measures the maximum distance 
of a point in a set X to the nearest point in Y or vice versa. 
Generalization to 3D uses mean or median over all slices. 
The measure indicates how much manual corrections are 
required. An ideal value of HD equal to zero reflects com- 
plete agreement of the contours. In addition, the number 
of slices with a HD greater than 3 mm is reported and 
compared to the total number of slices with prostate infor- 
mation. A threshold of 3 mm is often seen as clinically 
acceptable [22]. A second criteria is the volume overlap 
(or Dice similarity coefficient) defined by 



dyoiiX, Y) = 2 



\xnY\ 

\X\ + \Y\' 



(9) 



where | ■ | is the number of voxels contained in a region. 
Finally, accuracy is defined as 



d acc (X, Y) — I — 



\FP\ + \FN\ 
\TP\ + \FN\ 



(10) 



with TP = X Hi Y volume included in both X and Y (true 
positive), FN = X Hi (^Y) volume of X not included in 
Y (false negative) and FP = (^X) Hi Y volume of Y not 
included by X (false positive). Both values range from 0 
to 1, with optimal value 1. Volume overlap indicates how 
much of the prostate has been detected by the approach 
while accuracy shows how incapable the method is to 
select the true prostate pixels. 

In addition to the quantitative metrics, we have per- 
formed a small pilot test on 8 new patients comparing 
time expenditure using the proposed method and man- 
ual delineation. The time expenditure for the proposed 
method includes the definition of control points and the 
correction of the contour obtained by the method for 
each patient. The used mean shape model and variance 
was calculated from the training data as described above. 
Time measurements were obtained by two independent 
physicians for each case. Manual delineations, definition 
of control points and corrections were performed using 
the treatment planning system Eclipse™. 

Results and discussion 

The evaluation is performed on the deformed mean best 
fitting ellipses, i.e., on it(p \ CP) in formula (7). Additional 
evaluations are done in Additional file 1 for 7t(p \ S, CP). 
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Table 1 Evaluation metrics comparing BFE to manual delineations, and comparing deformed MBFE to manual 
delineations 



Test set 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 












BFE 












Dice 3D 


0.96 


0.97 


0.95 


0.96 


0.94 


0.94 


0.96 


0.96 


0.97 


0.93 


Accuracy 


0.93 


0.94 


0.91 


0.91 


0.88 


0.88 


0.92 


0.93 


0.94 


0.85 


HD mean 


1.32 


1.09 


2.24 


1.34 


1.94 


1.61 


1.41 


1.65 


1.59 


2.49 


#HD > 3mm 


0(6) 


0(8) 


3(12) 


0(8) 


2(11) 


0(13) 


0(9) 


0(15) 


1 (13) 


2(10) 












Deformed MBFE 












Dice 3D 


0.92 


0.93 


0.90 


0.91 


0.88 


0.84 


0.92 


0.88 


0.89 


0.88 


Accuracy 


0.84 


0.84 


0.81 


0.82 


0.74 


0.70 


0.84 


0.73 


0.80 


0.74 


HD mean 


2.38 


2.58 


6.12 


2.76 


3.79 


4.32 


2.73 


5.59 


4.88 


4.63 


#HD > 3mm 


1 (6) 


1 (8) 


9(12) 


4(8) 


7(11) 


8(13) 


2(9) 


13(15) 


8(13) 


9(10) 



(unit: Dice 3D and accuracy in percentage, HD mean in mm, number of slices with HD > 3mm compared to the total number of slices in brackets). 



Table 1 contains the distance metrics defined in (8) - (10) 
comparing the manual delineation and BFE for each test 
data set, and comparing the manual delineation and the 
deformed MBFE described by n(p | CP). The high Dice 
similarity coefficient and accuracy values and small Haus- 
dorff distances between manual delineations and BFE 
confirm the stacked ellipses model This is also reflected 
by the count of slices with a HD greater than 3 mm in each 
test data set. Only 7.6% of the in total 105 slices have a 
HD greater than 3 mm. The values show the best possi- 
ble description of the test cases by the proposed model. 
Furthermore, Table 1 presents the metrics comparing the 
manual delineation and the deformed MBFE for each test 
data set. The distance metrics reveal the fairly accurate 
results. The values indicate that the deformed MBFEs 
used as initial contours for final delineations will lower the 
time expenditure of the delineation procedure. However, 
more slices with a HD greater than 3 mm can be observed, 
59% of the 105 slices. Particularly, the test sets 3, 7 and 10 
have a ratio greater than 75%. Figures la to le illustrate 5 
slices of the BFE evaluation of test patient 3 from Table 1. 
Figures If to lj illustrate 5 slices of the deformed MBFE 
evaluation of test patient 3 from Table 1 with a volume 
overlap of the manual delineation line and the deformed 
mean shape of 0.90 and accuracy 0.81. 



Table 2 shows the median and median absolute devi- 
ation (MAD) for the data groups "test data" "training 
data" and "all data". A BFE volume overlap for all data of 
0.954 ± 0.010 (medianzbMAD) and accuracy of 0.908 ± 
0.020 confirm the model further (Table 2). Similar values 
in the subset of test data and training data are indicat- 
ing model robustness. In addition to the BFE results, 
Table 2 summarizes the results by median and MAD of the 
distances between manual delineation and the deformed 
MBFE. A median volume overlap of 0.899 ± 0.021 and 
accuracy of 0.807 ± 0.035 of the test data show further 
the power of the prior. The deformation of the prior is 
done by the control points and can be computed directly 
since there is no sampling or estimation involved at this 
point. 

To evaluate the robustness of the model, we randomly 
split 10-times the set of 33 patients into a training set with 
23 cases and a test set with 10 cases. Figure 4 shows the 
evaluation distances between the manual delineation and 
the deformed MBFE. The central mark is the median, the 
edges of the box are the 25th and 75th percentiles and 
the whiskers extend to the extreme data points. The figure 
shows only small variation between the different permu- 
tations, thereby demonstrating robustness of the stacked 
ellipses model. 



Table 2 Evaluation metrics comparing BFE/MBFE and manual delineations for different data groups 


Data group 






Test data 


Training data 




All data 


BFE 


MBFE 


BFE 


MBFE 


BFE 


MBFE 




median 


1.604 


4.052 


1.840 


3.806 


1.810 


3.806 


HD mean [mm] 


















MAD 


0.274 


1.305 


0.286 


0.427 


0.289 


0.538 




median 


0.959 


0.899 


0.954 


0.903 


0.954 


0.903 


Dice 3D [pet] 


















MAD 


0.008 


0.021 


0.007 


0.013 


0.010 


0.019 




median 


0.918 


0.807 


0.908 


0.800 


0.908 


0.800 


Accuracy [pet] 


















MAD 


0.015 


0.035 


0.014 


0.031 


0.020 


0.036 
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permutation permutation permutation 



(a) Hausdorff (b) Dice 3D (c) Accuracy 

Figure 4 Evaluation results between MBFE and manual delineations for 1 0 random permutations in training and test sets consisting of 23 
and 1 0 patients, respectively. The Hausdorff distance is depicted in (a), the percent-volume overlap measured by the Dice coefficient in (b) and 
the accuracy in (c). The central mark is the median, the edges of the box are the 25th and 75th percentiles and the whiskers extend to the extreme 
data points. 



Results from using MCMC to further optimize the 
delineation, as described in the Methods section, are only 
presented in Additional file 1 since a slight improvement 
comes at the cost of large computation time. 

The time comparisons indicated an average of 30% time 
saving using the proposed method compared to man- 
ual delineation. The time measurement of the proposed 
method includes the definition of controls points as well as 
the correction of the estimated contour by the physician. 

Conclusions 

The presented results demonstrate the potential of the 
proposed method in modeling the prostate by slicewise 
best fitting ellipses. Deformation of the mean shape using 
control points gives very good results with little compu- 
tational cost. Hence we believe that providing physicians 
with a good initial contour is beneficial in the clinical 
praxis of radiotherapy treatment. 

The corrections of generated delineations based on few 
control points were not streamlined in the workflow of 
the physicians, and the task of correcting contours is not 
part of their everyday activity. Furthermore, corrections 
were not done directly after the definition of the control 
points and sometimes by different physicians, and physi- 
cians had to deal with a different orientation of the data 
set in the treatment planning system than in the diag- 
nostic MRI. These issues must and can be solved for a 
well designed system. Therefore, a time saving of 30% 
likely represent a lower limit, and has to be validated 
in a well designed and properly powered study. Further- 
more, we expect larger time savings in data sets where 
the prostate is imaged in a higher number of slices. In the 
extreme case, if the prostate is visible in only three slices, 
the BFE approach would not give any benefit using the 
current control point method. The study of inter/intra- 
observer variability using the proposed method compared 
to manual delineation was considered to be beyond the 



scope of this study and is left open for interesting future 
work. 

In addition, the results show a precise description of 
the prostate by the BFE model with an average volume 
overlap of 95%. The high performance of the deformed 
mean shape model using the control points explains the 
small improvement by applying MCMC. Nevertheless, an 
improvement of the likelihood in the posterior distribu- 
tion or by an active appearance model [23] is a field of 
further research as elaborated in Additional filel. A clear 
disadvantage of an additional method like MCMC is the 
extra computation time. 

Further improvements can be achieved in the constraint 
and regularization terms, e.g., by considering the surface 
curvature versus changes of the ellipses parameters. We 
do not expect abrupt changes between neighboring slices 
around the central slice, but larger changes between slices 
towards the ends can be permitted, particularly in the 
length of the first and second principal axis. Also, the 
reduction of manual interaction in the proposed method 
is left for future work. 

Additional file 



Additional file 1 : Supplementary materials. Article containing i.) a 
detailed description of the relative coordinates systems (e.g. ICS, PCS) on 
the basis of the DICOM file structure, ii.) post-processing procedures as for 
example reordering and introduction of constraints, Hi.) a discussion of 
different control point method with construction of the transformation 
matrix A n dcp and the parameter rf l , iv.) elaboration of the posterior 
distribution, and v.) a section with additional data analysis. 
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